Massless particle creation in a f(R) accelerating universe 

S. H. PereiraJ*] J. C. Z. Aguilar, and E. C. Romao 
Universidade Federal de Itajubd, Campus Itabira 
Rua Irmd Ivone Drumond, 200 
35903-087 - Distrito Industrial II, 
Itabira, MG, Brazil 

^ ■ Abstract 



In this paper we present an alternative mechanism to generate dark matter particles through- 
out the decay of virtual scalar massless particles into real particle- antiparticle pairs. In order to 
illustrate such mechanism, we discuss quantitatively the problem of massless particle creation in a 



flat, homogeneous and isotropic universe expanding with a scale factor proportional to i 8 / 3 , which 
$— i , follows directly from a particular f(R) theory of gravity. We find that the spectrum of the total 

number and total energy of created particles can be determined exactly for any value of the phys- 
ical time in a matter dominated universe. We also find that the main contribution to the total 
■^J- . number of particles and total energy comes from small wavenumbers k and particles with large 

m 

values of k are produced only in the future. In the absent of an appropriated energy scale, if we 

OC 

choose the minimal mode k as corresponding to the maximum wavelength inside the Hubble hori- 
zon, we can estimate the present energy density of such massless particles as 10 3 times the present 
critical energy density. By applying the method of space-time correlation function we conclude 

X 

that such particles are virtual ones, forming a background field similar to the cosmic microwave 
background radiation. This opens the possibility that the decay of these virtual particles into 
massive particle-antiparticle pairs could be an alternative mechanism to generate particles of dark 
matter throughout the universe evolution. Moreover, this theory is naturally accelerated, giving 
an explanation for the current phase of acceleration of the universe. 
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I. INTRODUCTION 



The phenomenon of particle creation in an expanding universe has been studied by several 
authors lMldj after the pioneering work of Parker 11J. One of the most interesting results 
from Parker's work is that in a radiation dominated universe there is no creation of massless 
particles, either of zero or non-zero spin. In the description of particle production by the 
gravitational field, a widely used method is that of instantaneous Hamiltonian diagonaliza- 
tion [7| suggested by Grib and Mamayev [9|. Of particular interest are the works relating 
the gravitational particle production at the end of the inflation as a possible mechanism to 
produce super-heavy particles of dark matter 

All this searches are in the framework of standard general relativity. Recently, non- 
standard gravity theories have been proposed as an alternative to explain the present accel- 
erating stage of the universe only with cold dark matter, that is, with no appealing to the 
existence of dark energy. This is naturally obtained in f(R) gravity theories (for a review, see 
[3]). In such an approach, the curvature scalar R in the Einstein-Hilbert action is replaced 



by a 
case 



general function f(R), so that the Einstein field equation is recovered as a particular 
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22( | . A review of various modified gravities considered as gravitational alternative 



for dark energy, an unified description of the early-time inflation and the solutions of some 
related problems in f(R) theories have been presented by Nojiri and Odintsov 23(. In a 

I L 

recent work [24j] we have studied the quantum process of particle creation in a radiation 
dominated universe in view of a f(R) = R + (3R n theory, which has the correct limit of the 
standard general relativity for (3 = 0. In that paper we have showed that both massive and 
massless scalar particles can be produced by purely expanding effects, contrary to Parker 
results. 

In the present work we investigate the particle production in a f(R) = j3R n theory in a 
flat and matter dominated universe. We show that scalar massless particles can be produced 
as an possible alternative to create dark matter particles through the decay of these particles 
in particle-antiparticle pairs. Moreover, such f(R) naturally provide an accelerated stage 
for the evolution of the universe, without to appeal to a dark energy component. 

Although a scaling factor of type t 8//3 seems to be artificial, this model allows an exact 
solution for the spectrum of created particles in terms of the real physical time t and with 
initial conditions at t{. Thus, the treatment presented here can be the basis for future 
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research related to the creation of matter in the universe for various models of evolution. It 
is important to note that most studies in the literature consider the creation process only 
in limiting cases, t — > oo. Here we presents finite results concerning real and finite time t. 



II. GENERAL THEORY OF SCALAR PARTICLE CREATION 

The canonical quantization of a real minimally coupled scalar field in curved backgrounds 
follows in straight analogy with the quantization in a flat Minkowski background. The grav- 
itational metric is treated as a classical external field which is generally non-homogeneous 
and non-stationary. The basic equation for the study of scalar particle creation in a spatially 
flat Friedmann- Robertson- Walker geometry is (J): 

Xfcfo) + <4(v)xk{v) = o , (i) 

with 



"l(v) = k2 + m lff and m lff =m 2 a 2 (ri) - , (2) 



where k is the Fourier mode or wavenumber of the particle, Uk and m e ff represents the 
frequency and the effective mass of the particle, respectively. a{rj) is the cosmological scale 
factor in terms of the conformal time rj and the prime denotes derivatives with respect to 
it. The conformal time rj and the physical time t are related by 

Following standard lines, the quantization can be carried out by imposing equal-time 
commutation relations for the scalar field \ and its canonically conjugate momentum 7r = 
X', namely [x(x, fj) , ir(y, rj)] = i8(x — y), and by implementing secondary quantization. 
After convenient Bogoliubov transformations, one obtains the transition amplitudes for the 
vacuum state and the associated spectrum of the produced particles in a non-stationary 
background [3|, |4j • 

Usually, the calculations of particle production deal with comparing the particle number 
at asymptotically early and late times, or with respect to the vacuum states defined in two 
different frames and do not involve any loop calculation. However, the main problem one 
encounters when treating quantization in expanding backgrounds concerns the interpretation 
of the field theory in terms of particles. The absence of Poincare group symmetry in curved 



space-time leads to the problem of the definition of particles and vacuum states. The problem 
may be solved by using the method of the diagonalization of instantaneous Hamiltonian by a 
Bogoliubov transformation, which leads to finite results for the number of created particles. 

To proceed, note that Eq. ([!]) is a second order differential equation with two independent 
solutions. Each solution Xk must be normalized for all times according to 

W h {rj) = Xk(v)x*' k (v) ~ X'MX*M = -2i , (4) 
and also they must satisfy the initial conditions at the time r\i. 



XkiVi) = l/V^M , x'kivi) = W^iVi) • ( 5 ) 

The Bogoliubov coefficients can be calculated and a straightforward calculation leads to 
;he final expression for the total number of created particles and antiparticles in the k mode 

N k ( V ) = N k (v) = ^j\Xk(v)\ 2 + ^\x k (v)\ 2 ~\- (6) 

Another important quantity is the total energy per mode, E k {r\). It is proportional to the 
product of the angular frequency u>k with the total number of created particle in that mode, 

E k {rj) = 2u k (r])N k (r)) (7) 

where the factor 2 stands for the fact that the energy of particles and antiparticles are equal 
in this case. 

Thus, the total number density of created particles n and the total energy density e are 



readily obtained by integrating over all the modes 



6] 



1 f°° 

i r°° 

All the expressions ([6])- (Q can be rewritten in terms of the physical time t by using the 
relation (131). 



III. PARTICLE CREATION IN A SIMPLE f(R) THEORY 



Now we will apply the above results to the study of massless scalar particle creation in a 
matter dominated and accelerating universe with scale factor of the form a(t) ex t 8 / 3 , which 
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follows readily from a simple f(R) theory, discussed in detail by Allemandi et al 16], namely 



f(R) = m n . (io) 

Taking into account the dominant energy condition, R n must be positive definite for even 
integer values of n and /3 should be fixed as (3 > for n < 2 and < for n > 2. 

Now let us consider the case of matter dominated universe (oj = 0). The modified 

n 

Friedmann equations can be obtained [16j (apart from integration constants) 



a{t) 



3e ^ n/s 



2n(3 - n) 



SnGpo 



1/3 / £ \ 2n/3 



to/ 



(11) 



where e = 1 for odd values of n and e ± 1 for even values. The dot denotes derivatives with 
respect to physical time t and p is the present day value of the energy density. Note that 
the case n = 3 is singular. With this expression we can calculate the deceleration parameter, 

We see that accelerated solutions (g < 0) are obtained in the cases ri<0orn>3/2. 

In order to illustrate the phenomenon of particle creation related to this theory we must 
to set particular values to the parameters n, (3 and e. We choose n — 4, which corresponds 
to an accelerated a(t) according to (fT2|) . (3 < and e = — 1. With such parameters, the scale 
factor is 

and the deceleration parameter q = —5/8. The Hubble parameter H = a /a = 8/3t deter- 
mines the actual age of the universe as to = 8/3Hq — 18.3 x 10 9 years. It is interesting to 
note that it does not depends on the (3 parameter. 

To follow up we need the scale factor a(t) in terms of the conformal time i]. This can be 
obtained by integrating ()3]) and rewritten a(t) in terms of rj, 

/ \ 8/5 ^ 

a (^) =a of^j ^75? 0<r ? <oo. (14) 

In terms of the conformal time, the early universe (the past) corresponds to rj — > oo, and 
the late universe (the future) corresponds to rj — > 0. 

Before we proceed with the calculation of the spectrum of non-massive particles created 
during the evolution of the universe, we must choose a value for the /3 parameter. We will 



choose (3 such that a(t ) = oo = 1, for t = 8/3H , i.e., the scale factor is normalized to 
unity for the present time. To satisfy this condition we must have 



3 X12 



/? = -4ttG - p . (15) 



This gives /3 ~ —6.6 x 10 86 GeV 2 if we set po of the order of the critical density. Thus the 
scale factor is 

/ o \8/5 -, 

"to-isnd ^' (16) 

It is important to note that we can recover all the expressions in terms of the physical time 
t by the substitution 



3to ft_ x ~ 5/3 



'=ssUJ ' (17) 

The equation (j2J) for the mode function in the case of a massless particle (m = 0) is 

104 



k z - 



Xk(v) = . (18) 



25r/ 2 _ 

In order to satisfy the initial conditions (jSJ), the solution of this equation is given by 

Xk (v) = A(k, Vi)VvYzL (kv) - B(k, riVvty (kv) , (19) 

where J u and Y u are the Bessel functions of first and second kind, respectively, v is its order 
and A and B are complex constants depending on k and on the initial time T]i, 

8Jai(%) - 13 J ii (kr]A + 1 21J 21 (kr/i) 



A(k,r]i 



21v^[j2i(kr)i)Yu(kr) i ) - Yk (krji) Ju (%)] 



8F 3 i (kr)i) - 13Yn (km) + 1 2\Y%x (km) 

B(k,m)= 15 75 (20) 

21v^[J^(%)Yii(A;77i) - Y2i(%)Ju(%)] ' 

10 7 10 10 v ' ' 10 J 

where / = a/— 1. 

By fixing m and taking some particular values of k and 77 it is easy to check that this 
solution correctly satisfy the normalization condition (jl]). 

The solution presented in ( fl9|) is exact and, at first, it may be used to calculate the total 
number density of particles OH]) and the total energy density for any instant of time rj 
after the initial time m- To recover the expressions as a function of physical time t we use ([3]) 
and all quantities can be described in terms of t. In fact, the final analytical expression for 
Nk is quite complicated, but can be obtained in relatively simple manner with the help of an 
algebraic manipulation software. Thus, the exact behavior and several interesting features 
of the functions ([6]) - © can be obtained by graphical analysis of the spectrum of created 
particles. This is what we will do in the next section. 
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FIG. 1: (a) - Spectrum of the total number and (b) - the total energy of created particles as 
function of the physical time t in units of Hq 1 for some different values of the mode k. The present 
time is indicated by the vertical dashed line, to- 

IV. SPECTRUM OF CREATED PARTICLES 

Now let us study the evolution in time of the total particle number N/~ and the total 
energy per mode Ek, and try to estimate n and e for such particles. 

First we need to establish an appropriate time scale in order to make the analytical study. 
As we have seen, the present time corresponds to t = 8/3H = 5.8 x 10 17 s, so that an 
appropriate time scale is in terms of -f^o" 1 - Using ( I17J) we can estimate this time as i]q = 1.6. 
Now we will establish the value of the initial time 7/j. As we are dealing with a universe 
dominated by matter and we know 26| that the matter and radiation decouple at about 
ti = 10 13 s ~ 4.5 x lO -5 -/?^ 1 , by using (fiTl) we can estimate r]i = 1.4 x 10 8 (remember that, 
in the conformal time, the absolute past corresponds to rj — >■ oo). Thus we are interested in 
the interval 1.6 < 77 < 1.4 x 10 s , which corresponds to 4.5 x KT 5 //^ 1 < t < (S/S)//^ 1 . 

Now we can replace (fl9|) into (J6]) to obtain the spectrum of the total number of particles 
created for each mode k. With this we use ( TP?]) to write the result as a function of physical 
time t. The spectrum of the total number of created particles for some values of the k 
mode is represented in the Figure l-(a). By analyzing the Figure, we can see that for every 
mode k there is a critical value of the time t for which the number Nk of created particles 
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grows abruptly and diverges. After that the creation process stops for the corresponding 
mode. Mathematically this occurs because in the expression fl6]) the frequency Uk appears 
in the denominator. Let us express in terms of the physical time: 



We can see that when k 2 = 3 1 / 3 (351/8192)t 10 / 3 we have oo k = and the first term of flUJ) 
diverges. Physically, the significance of this divergence is that, for a given mode k, the 
values of t > (8192/351) 3 / 10 3~ 1 / 10 £; 3//5 turns the frequency squared negative, io\ < 0, and 
consequently the state of minimum energy and the quantum vacuum are not well-defined and 
the creation process ceases for these modes. Considering the present time in units of Hq 1 , 
the limit of k is k ~ 1.2747549. This indicates that up to the present time particles with 
modes k < k does not exist. We will see later that such modes corresponds to wavelength 
greater than the Hubble horizon. The modes with k > ko started to be created in the past 
and the divergence will occur only in the future. 

Note that an infinite number of particles per mode can be created until the process to 
stop. In the Figure 1, for instance, it is represented by the modes k = 0.5 and k = 1.0, 
where the divergence occurs in the past. For the mode k = 1.27 the divergence is occurring 
today. For the modes k = 2.0, k = 3.0 and k = 10 a significance amount of particles will be 
created only in the future. We see that the greater the value of k, the divergence will occur 
in the future, meaning that particles of higher energies can be created only in the future. 

Is it really divergent the total number of created particles per mode? Mathematically 
the expression ([6]) tells us that its happens due to the factor u k in the denominator. But 
perhaps there is some mechanism capable of interrupting the process of creation before it 
goes to infinity, so that the total number of created particles per mode is finite. We believe 
that such hypothesis is true based on the spectrum analysis of the total energy of created 
particles per mode. The expression ((7j) contain the factor Uk multiplying the total number 
of particles Nk, just canceling the factor Uk in the denominator of the first term of Nk- So 
we can conclude that the total energy of created particles per mode k must be finite. In fact, 
the graphical analysis shows that. In Figure l-(b) we present the energy spectrum of created 
particles per mode as a function of the physical time t. We clearly see that the energy grows 
to a maximum value for each mode k but it does not diverges. For larger values of the time, 
the frequency becomes negative and the creation process stops, as discussed earlier. 




(21) 
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FIG. 2: (a) - Spectrum of the total number and (b) - the total energy of created particles as 
function of the mode k for some different values of the time, where to represents the present time. 

Another interesting analysis that we can do based on the analytical results is the study 
of the behavior of Nf~ and Ef~ as a function of k for some fixed time. This is shown in 
Figure 2 for three different instants of time, namely t = to/2 (past), t = to (present) and 
t = 2t (future). For N k , as expected, the total number of particles per mode appears to 
diverges for a given value of the mode k, but decreases with increasing k. The same goes 
for E k , although the total energy per mode is always finite, as expected. The decrease of 
these values with the increase of the mode k is an interesting feature, indicating that higher 
modes contribute much less to the total number of particles as for the total energy. 

We have seen that the total number of particles is infinite, and the same will happen with 
the total number density of particles ([8]), unless some mechanism stop the creation process 
before the divergence occurs. We believe that this hypothesis is true, since the total energy 
per mode is finite. Let us estimate the total energy density ([9]) for the present time to- Since 
we do not have a mass scale, we need an appropriate energy scale for the modes in order to 
do the integration. 

From the above analysis we conclude that for the present time the value of the lowest 
mode is k ~ 1.2747549. For lower values of k the frequency becomes negative and there 
is no creation of these modes. We have also seen that the contribution to the total energy 
decreases with increasing k, so the integration in ([9]) can be truncated at some value of k in 
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the case of a numerical integration. First remember that the mode k is related to physical 
wavenumber k p h by k = a(t)k p h. The physical mode k p h is related to the physical wavelength 
^ph by kph = 27[/X p f l . Therefore a minimum value of k corresponds to a maximum value of 
the physical wavelength X p h, 

k mm = a (t)- , (22) 

\ ' \max ' \ ' 

We choose the maximum value of A as the size of the Hubble horizon today, \™ h ax = cto = 
c8/3Hq ~ 3.52 x 10 28 cm. With this we are requiring that only the modes with wavelengths 
shorter than the actual size of the universe contributes to the total energy, which seems 
plausible, since we are restricting ourselves to the modes inside the Hubble horizon. Modes 
with longer wavelengths have k < ko and are not allowed. Thus we have an appropriate 
scale of energy to estimate the total energy density of the particles being created at present. 
The minimum value of the physical wavenumber is k™/} 71 = 2ir/ct ~ 3.5 x 10 _42 GeV. With 
this we make a change of variables in ([9]) such that k$ corresponds to k mm and the energy 
density can be calculated numerically for the present time as 

1 f°° 

£o = tt^tt / k 2 E k (t )dk « 1.91 x 10" 44 GeV 4 . (23) 

This result is very interesting, since the value of the energy density of the particles being 
created at the present time corresponds to 10 3 times the value of the critical energy density, 
Po = 4.0 x 10~ 47 GeV 4 , forming an energetic background field similar to the cosmic microwave 
background radiation field. So, if we consider that part of these non-massive particles that 
are being created today are decaying into massive particle-antiparticle pairs, they may have 
enough energy to give rise to the dark matter in the universe. 

To finish, let us apply the method of a space-time correlation function suggested in 
in order to distinguish if such non-massive created particles are real or virtual ones. This 
method was generalized by Pavlov |7| for a scalar field with nonconformal coupling. It is 
easy to check that for the present time the metric (1141) is changing adiabatically for modes 



k » ko, 



k 



«1, (24) 

and furthermore we have a'(r]i) — > for r/i — > oo, so that all the analysis of Q are valid 
here, and we conclude that such a pairs of created particles should be interpreted as virtual 
particles. As they are virtual particles, its total energy does not contribute to the total 
energy of the universe. 
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V. CONCLUDING REMARKS 



We have addressed the problem of minimally coupled massless particle creation in a 
f(R) = PR 4 theory for a matter dominated universe. We have verified that this particular 
type of f(R) correctly reproduces an universe that evolves in accelerated manner (q = —5/8). 
In addition we have showed that an infinite amount of non-massive virtual particles with 
mode k < k has been created in this particular model, although we believe that some 
mechanism stops the creation process before the number of total particles goes to infinity. 
Such an idea is based on the fact that the total energy per mode is finite, so we think that the 
total number of particles is also finite. Making an appropriate choice of the energy scale for 
the minimum mode k as those corresponding to the possible wavelengths within the Hubble 
horizon, the estimate for the total energy density of the virtual particles pairs being created 
today shows that the value is about 10 3 times larger than the critical energy density. As they 
are non-massive virtual particles, its total energy does not contribute to the total energy of 
the universe, but on the framework of the quantum mechanics, these particles could decay 
into massive real particle-antiparticle pairs, as occur in pair production of electron-positron 
in the decay of an energetic photon. This opens the possibility that these virtual particles 
being created today could give rise to the dark matter in the universe. In order to be 
a source of dark matter particles, the mechanism of decay of these particles into massive 
particle-antiparticle pairs needs to be studied more carefully. 
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